Complexity of 2D random laser modes at the transition 
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The spatial extension and complexity of the eigenfunctions of an open finite- size two-dimensional (2D) ran- 
dom system are systematically studied for a random collection of systems ranging from weakly scattering to 
localized. The eigenfunctions are obtained by introducing gain in the medium and pumping just above thresh- 
old. All lasing modes are found to correspond to quasimodes of the passive system, for all regimes of propaga- 
tion. We demonstrate the existence of multipeaked quasimodes or necklace states in 2D at the transition from 
localized to diffusive, resulting from the coupling of localized states. 

PACS numbers: 42.25.Dd, 42.55.Zz 



Transport in random media is driven by the nature of the 
underlying eigenmodes. Propagation is diffusive when the 
modes extend spatially, while spectral level overlap occurs in 
transmission spectra. As the degree of the overlap decreases, 
transport is inhibited and modes become spatially localized 
ill]. The theory of Anderson localization predicts a transition 
between localized and extended eigenstates for spatial dimen- 
sions larger than two [2]. Renewed interest in the localiza- 
tion transition has been boosted by the active ongoing search 
for localization of Bose-Einstein condensate in laser speckle 
fields H,^, 5], the recent observations of the slowing down of 
diffusion in ultrasounds [6], microwave itzt] and time-resolved 
optical Isl] experiments, and new theoretical progresses 
towards an analytical description of the metal insulator tran- 
sition (MIT). The question of the spatial extent of the modes 
near the Anderson transition is also central in random lasers 
where the threshold is still high. The threshold may 
vary by orders of magnitude between localized systems where 
the modes are spatially confined and diffusive systems where 
the modes are extended. 

If the non-coexistence of extended and localized states was 
recently confirmed numerically [13], recent experimental ob- 
servations of necklace states in nominally localized optical 
and microwave iflsll one-dimensional layered systems 
have attracted much attention. Modes overlapping both in 
space and frequency may couple and form multipeaked ex- 
tended states even in the localized regime lITdl . Although 
scarce, they are predicted to play an outsized role in trans- 
port fivll in contrast to isolated localized states |fl8|]. Pendry 
argued that the picture of ID necklace states should generalize 
to 2D and 3D localized random media fl7]. However, besides 
earlier calculations in a percolation model, no observations of 
necklace states in higher dimensions than one were reported 
|[l9[ |20^ . In this context, the marginal dimension of two is of 
special interest since a pseudo-MIT is possible in 2D finite- 
sized opened samples where transport can be ballistic, diffu- 
sive (in contrast to ID) or localized [21]. Besides, the spatial 
profile of the modes is readily accessible as opposed to 3D, 
as shown in a recent observation of localized states within a 
2D random microwave cavity ilTA . The evidence of necklace 



states would therefore shine a new light on the Anderson tran- 
sition, suggesting a possible scenario at the transition where 
coupling between increasing number of states could lead to 
a quantum percolation-type of transition [23]. This is to be 
compared with the filament-like fractal picture of the modes 
at the transition suggested by Aoki [24], which lead to the ex- 
tensive use of multifractal analysis to explore the MIT 1I25I1 . 

Fundamental changes in the nature of the quasimodes are 
also expected in open random systems, as their spatial exten- 
sion grows beyond the sample dimensions and their linewidth 
broadens with increasing leakage. This is analogous to quan- 
tum systems with an increasing degree of opening where 
the wavefunction becomes complex- valued with its standing 
wave component being progressively replaced by a compo- 
nent traveling toward the opened boundaries. The intensity 
distribution of the spatial pattern of the mode is expected 
to cross over from Porter-Thomas to Rayleigh distribution 
when new leakage channels ope n ll27ll . The complexity fac- 
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together with the phase 
rigidity, p = {1 — cf') / {1 ^ q^) ll29l], was introduced in the 
field of quantum chaos to quantify the degree of complexity 
of the eigenfunction, ^, and the mutual influence of neigh- 
boring resonances 11301] . The complexity factor crosses over 
from for real standing wavefunctions to 1 for purely travel- 
ing waves. To the best of our knowledge, it has never been 
used as a probe to characterize the transition from localized to 
extended states in finite disordered systems. 

In this letter, we explore the nature of the quasimodes of 
2D random media when scattering strength is increased. To 
obtain the quasimodes, the passive random system is embed- 
ded in a gain matrix and pumped just above the threshold. 
Single-mode lasing with resonant feedback is observed in any 
regime between extended and localized, with lasing modes 
corresponding systematically to a quasimode of the passive 
system. The spatial extension of the quasimodes, their com- 
plexity factor and the lasing threshold are calculated for a sta- 
tistical ensemble of random configurations for each value of 
the scattering strength. These quantities reveal the change of 
regime and the transition from diffusive to localized. A de- 
tailed analysis of the phase probability distribution for each 
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mode reveals multipeaked wavefunctions in the vicinity of the 
transition when the locaHzation length is comparable to the 
sample size. These extended multipeaked quasimodes corre- 
spond to coupling of localized isolated states. 

We consider a two-dimensional random collection of paral- 
lel dielectric cylinders with infinite extension, radius r = 60 
nm and refractive index n, embedded in a background ma- 
trix of index 1. Volume fraction is = 40 % and system 
size is = 5x5 /im^. Maxwell equations for TM polar- 
ization are modeled using the finite-difference time-domain 
method [31]. Open boundary conditions are approximated by 
perfectly matched layer absorbing boundaries. The index of 
refraction n is varied from 1.05 to 2.0, in step of 0.05, cor- 
responding to scattering mean free path ranging from 50 /im 
to 0.1 /im. For most of this range, modes are short lived with 
strong spectral overlap, preventing individual excitation of a 
mode at its eigenfrequency by a monochromatic source in or- 
der to obtain its wavef unction. We showed recently that, when 
operating just above threshold, the first lasing mode of an ac- 
tive random system corresponds to a quasimode of the passive 
system, even in weakly scattering systems ll32|] . Introducing 
gain and pumping just above the threshold is therefore an al- 
ternative to select a quasimode of the passive cavity and obtain 
its spatial distribution, even in a system where modal overlap 
dominates. To model the gain, we couple the population equa- 
tions of a four level atomic system to the Maxwell equations 
via the polarization equation 1I33I1 . The gain naturally selects 
the mode with the longest lifetime and the best spectral over- 
lap with the gain curve. 

We study 150 random configurations, with an average of 
10 per value of refractive index, n. For each configuration, 
the pump power is adjusted just above the threshold to ob- 
serve single-mode lasing. Following the method using the 
spatial correlation introduced in [32], we check that the las- 
ing modes are the quasimodes of the passive system which 
provide the necessary feedback for lasing. Six spatial distri- 
butions of quasimodes corresponding to increasing values of 
the refractive index, n, are shown in Fig.[TJ illustrating differ- 
ent degrees of spatial extension of the wavefunctions within 
the system. As modes extend toward the system boundaries 
with decreasing refractive index, they acquire an imaginary 
part corresponding to a traveling wave component, in contrast 
with localized modes, which are purely standing waves. This 
is seen in the phase probability distribution between and 27r. 
As shown in insets in Fig.[T] the phase probability distribution 
is peaked around and tt when the mode is localized (n=2), 
while it is more uniformly distributed in the extended case. 
Note that for values of the refractive index as low as 1 .05, scat- 
tering is weak and the field is rather concentrated at the edges 
of the system. In that case, residual reflection either at the 
boundaries or from the PML layers may not be negligible and 
result in periodic patterns, similar to those of a Fabry-Perrot 
cavity, as seen in the first frame of Fig. [T] We checked that 
above n=1.10, this effect is insignificant and lasing is solely 
due to multiple scattering within the system. 

The two main features observed as the scattering strength 
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FIG. 1 : Spatial distribution of the magnitude of the lasing mode just 
above the threshold, as the refractive index of the scatterers is in- 
creased from n=L05 to 2.00 (random configurations are not neces- 
sarily identical). Each frame shows in inset the phase distribution 
between and 27r of the corresponding quasimodes. Note the double 
peaked distribution for n=1.85. 



is reduced are the spatial expansion of the eigenfunctions and 
the increase of their imaginary parts associated with leakage 
at the open boundaries. To quantify these two characteristics, 
the quasimodes are described in terms of their spreading fac- 
tor and their complexity factor. We define the spreading factor 
as ?7 = J J A{r)A(r')\r — r'\^dPrd?r' , where the field am- 
plitude A{r) is normalized A (r) = A{r)/[j J A'^{f)£r\^^'^ . 
The spreading factor is itself normalized to unity for uniform 
spatial distribution of the field. It measures the degree of spa- 
tial extension of the energy within the system alike the par- 
ticipation ratio for instance. However, due to the weighting 
factor |r — r'p, it also enables us to distinguish systems with 
spatial localization of energy inside the system, as in Fig. [T] 
for n=2, from concentration of energy at the boundaries of the 
system, as in Fig. [T] for n=1.05. Indeed, it can be less than 
1 for spatially localized modes, or larger than one when en- 
ergy is distributed near the system edges. This is reminiscent 
of distributed feedback lasers 1 34] where in the over-coupled 
regime (corresponding to 7^ < 1), energy is concentrated in- 
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FIG. 2: Color online, (upper) Complexity factor, and (lower) 
spreading factor, 77, averaged over sample configuration versus scat- 
tering mean free path, I. The fluctuations around the average, 
=t(^^ — (9^))"^^^^ and ±(?7^ — (77))"^/^, are represented by the bars. 
Inset: Localization length ^ versus index of refraction, n, calculated 
from the averaged lasing threshold (full line) and from independent 
scattering theory (full line where relevant (^ < L/2), dotted other- 
wise). The dashed line represents the mean free path, ^. The hori- 
zontal dotted line correspond to ^ = L/2 



side the laser as a result of strong feedback from scatterers, 
while in the under-coupled regime (7/ > 1), energy is con- 
centrated at the edges since the lasing modes result from scat- 
tering at the boundaries in order to maximize the gain vol- 
ume 1 35]. The spreading factor and the complexity factor 
are computed for each mode [36] and averaged over sam- 
ple configurations for each value of n. They are shown in 
Fig. [2] as a function of the scattering mean free path, ^, cal- 
culated using Mie theory for infinite cylinders of refractive 
index n. As the complexity factor explores values between 
0.78 and 0.06, it shows clearly two different regimes, with 
a crossover around ^=0.14 jim corresponding to n=1.8. A 
transition occurs around the same value of I for the spread- 
ing factor, which ranges between 1.23 and 0.16. This value 
certainly is not universal and depends on the sample size but 
it should correspond to a localization length, ^, of the order of 
the sample size. We confirm this hypothesis by calculating ^ 
directly from the value of the lasing threshold averaged over 
sample realizations, (P). Indeed, the spectral width of the 
modes, F, resulting from leakage at the boundaries is given by 
F = Fo-exp(— L/^) [22]. It is also directly proportional to the 
lasing threshold, P, since at threshold, losses are compensated 
by gain. Therefore ^ = L/(ln(Po) — ln(P)),where (P) des- 
ignates the average over sample configurations for each value 
of n and In (Pq) =20. 15 is obtained by extrapolating (P) at 




FIG. 3: Phase spatial distributions for phase values around (a) 0.947r 
and (b) I.OGtt. 

n = 1. The dependence of ^ on refractive index, n, is shown 
in the inset in Fig [2] and is compared to the theoretical expres- 
sion in the limit of independent scattering, as given in [ 22] by 
6^ = ^exp[7ri?e(/ce//)^/2], where ^ is the mean free path 
and kef / is the effective wavenumber. Both curves join at the 
crossover, n=1.8, which corresponds to ^ ~ L/2 =2.5 /im, 
where the two expressions for ^ and £^th start to be valid. 

Also shown in Fig [2] are the fluctuations of the com- 
plexity factor, — (^^))^^^, and the spreading factor, 
±(77^ — (7?))^^^. A significant increase of these fluctuations 
is seen at the crossover. We find that these larger fluctuations 
are correlated with the occurrence of peculiar phase distribu- 
tions around n=1.8, such as the one shown in the inset of 
Fig. [T] for n=1.85. Two narrow peaks are distinctly seen in 
the phase distribution at 0.947r and I.OGtt. The spatial distri- 
butions of the phase for values comprised in a tt/IO window 
around each of these two peaks are shown in Fig. [3^ and b. 
These two distributions delimitate two distinct spatial regions 
associated with standing components of the mode, which os- 
cillate at the eigenfrequency of the mode but with a phase lag, 
5(j) = 0.127r. This phase lag suggests that this mode results 
from the coupling between two distinct modes localized on 
each of the regions of Fig. [51 This would be the analog of the 
symmetric or anti- symmetric solutions to the coupled oscil- 
lators problem in the presence of leakage, which introduces a 
phase lag different from or tt between the components of the 
hybridized mode. To identify each of the two components of 
the original double-peaked mode of Fig. [T](n= 1.85), we use a 
perturbation approach where scatterers are removed in one of 
the two regions of Fig. [3] to selectively separate the two con- 
tributions. Each perturbed system is excited at the resonant 
frequency of the original unperturbed mode. The correspond- 
ing field distributions are shown in Fig.|4^ and b. They repro- 
duce the local features of each peak of the mode of Fig.[TJ but 
extend far beyond. The corresponding phase distributions are 
now single-peaked. Note that layers of randomly distributed 
scatterers (not shown) were added at the boundaries of the per- 
turbed system in order to increase the lifetime of the mode of 
Fig.H^, which would be impossible to excite otherwise due to 
strong leakage at the boundaries. 

We find that hybridized quasimodes with two or three peaks 
arise almost systematically at n=1.8 and 1.85. The modes ob- 
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FIG. 4: Spatial distribution of the magnitude of the quasimode to- 
gether with its phase distribution (in inset) for two different local 
perturbations of the random system of Fig. [T] for n=1.85, when ex- 
cited at the eigenfrequency of the original mode. The locations of the 
removed scatterers are shown by the circles. 

served here demonstrate the existence of necklace states in 
2D random media. They may play a crucial role at the tran- 
sition from localized to diffusive. Two different mechanisms 
were proposed, which may coexist |0]. One possible mech- 
anism is a gradual spatial expansion of the mode as scattering 
strength is diminished, with a progressive increase of the lo- 
calization length. In a different scenario analogous to a per- 
colation process, coupling of localized states could lead to 
extended structure as the spectral density of states increases 
to form necklace states The multipeaked phase distri- 

butions and correspondingly multipeaked modes found here 
with high probability around the transition, suggest that perco- 
lation is a possible scenario leading to a rapid crossover from 
localized to diffusive regimes in random media. 
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